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ABSTRACT 


Analytic  models  of  communication  networks,  combining  the  desired 
degree  of  simplicity  and  accuracy,  are  noticeable  by  their  absence 
from  the  literature.   Modern  networks  are  so  large  that  the  combina- 
torial possibilities  impose  an  overwhelming  burden  on  mathematical 
model  formulation.   High  speed  digital  simulations  provide  one  of  the 
only  means  of  solution  for  the  systems  planner. 

Military  communication  networks  offer  more  hope  for  analytical 
solution.   Much  smaller  than  their  civilian  counterparts,  most  military 
systems  can  be  modelled  satisfactorily  as  aggregations  of  queues. 
This  paper  provides  approximations  to  the  probability  of  blocking  in 
2  terminal,  k-channel  systems.   The  probability  of  blocking  is  defined 
as  the  probability  that  a  vacant  channel  is  available  for  connecting 
a  caller  to  his  desired  party  whenever  the  call  is  placed. 
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CHAPTER   I 
FORMULATION  OF  THE  PROBLEM 

1.   Introduction. 

Communication  systems  represent  an  integral  part  of  progress  in 
the  social,  business  and  technical  fields.   Certainly  this  statement 
needs  no  qualification  when  applied  to  the  military.   With  the  increase 
in  destructive  power  and  improvement  of  accuracy  in  modern  weapons, 
large-scale  dispersion  of  troops  on  the  battle  field  has  become  part 
of  the  doctrine  of  warfare.   Wider  separation  has  imposed  difficulty 
in  controlling  and  directing  forces   at  all  levels  of  command  from 
squad  leader  to  army  commander.   Rapid,  secure  and  reliable  communica- 
tion is  indispensable  in  marshalling  resources  to  obtain  the  military 
and  political  objective.   Indeed,  quality  and  quantity  of  communica- 
tions systems  often  provides  the  margin  of  victory. 

Military  communication  systems  encompass  a  wide  range  of  sophis- 
tication; the  spectrum  runs  from  the  hand-held,  short  range,  single 
channel  radio  to  the  most  modern  orbiting  satellite  system.   By  far 
the  most  important  in  terms  of  volume  of  traffic  and  ease  of  control 
is  the  multi-channel,  wide-band  radio  system  in  the  Army  division. 
This  system  closely  parallels  the  long-haul,  inter  city  trunking 
system  operated  nationally  by  the  Bell  System  and  independent  telephone 
companies.   Both  of  these  systems  have  a  network  of  radio  links  inter- 
connecting nodal  or  terminal  points. 

In  the  Army  division  system,  subscribers  are  tied  by  2-wire  cables 
to  the  switching  center  at  the  terminal.   Local  calls  are  switched  and 


connected  to  the  called  party  at  the  local  terminal.   Calls  to  other 
terminals  are  transmitted  over  radio  links  to  the  desired  terminal 
and  are  then  connected  to  the  2-wire  cable  of  the  desired  party. 

Civilian  networks  are  operated  with  automatic  switching  equipment 
to  route  calls  while  military  systems  are  at  present  using  manual 
methods  in  tactical  operations.   Another  difference  between  the  two 
systems  is  in  the  distribution  of  calls.   The  civilian  system  has  the 
bulk  of  calls  contained  within  the  local  area  whereas  the  military 
network  has  the  vast  majority  of  calls  going  outside  the  local  terminal. 
This  is  because  subscribers  in  the  military  system  are  located  in  close 
proximity  to  the  other  subscribers  in  the  same  terminal. 

Much  has  been  done  on  various  aspects  of  communication  networks. 
Perhaps  the  foremost  source  of  analysis  associated  with  developments  ■ 
in  this  field  is  the  Bell  System  Technical  Journals.   In  these  analyses 
efforts  have  been  made  to  arrive  at  mathematical  models  admitting  an 
accurate  display  of  the  significant  features  of  communications  systems. 
Practical  results  have  been  disappointing.   The  sheer  size  of  a  civil- 
ian system  is  staggering.   Millions  of  subscribers  tied  into  terminals 
that  are  connected  by  radio  links  of  up  to  a  thousand  channels  create 
immense  combinatorial  problems.   Consequently,  most  work  of  a  practical 
nature  in  the  analysis  of  the  operation  of  networks  is  done  by  computer 
simulation.   Ease  of  solution  provides  ample  compensation  for  the 
inelegance  and  inefficiency  of  the  method.   Military  systems  are  more 
amenable  to  analytic  solution  because,  as  will  be  seen  in  Section  1-3, 
radio  links  usually  have  a  maximum  of  24  channels  and  rarely  does  a 
terminal  interconnect  with  more  than  two  others.   Accordingly,  this 


paper  will  attempt  to  analytically  model  the  physical  operation  of  a 
military  communications  system. 

A  useful  parameter  of  a  communications  net  is  the  probability  of 
blocking,  P,  .   This  is  defined  as  the  probability  that  a  call  cannot  be 
routed  so  that  the  desired  party  is  reached,  because  all  available 
channels  of  communication  are  in  use.   The  theory  of  Markov  processes 
and  queues  was  used  to  obtain  the  value  of  this  parameter. 

Approximate  solutions  for  P,  will  be  obtained  for  a  2  terminal, 
k-  channel  system.   In  addition,  an  intuitive  approach  to  the  solution 
of  a  3  terminal  system,  again  interconnected  with  k-channel  links, 
will  be  given.   The  solution  for  the  2  terminal  model  will  be  preceded 
by  a  thorough  investigation  of  a  basic  system. 

2.   Organization. 

As  mentioned  above,  this  paper  will  attempt  to  develop  an  analytic 
model  for  finding  blocking  probabilities  for  small  communications 
systems.   An  orderly  progression  has  been  attempted  by  initially  treat- 
ing a  very  basic  model  and  then,  by  adding  real  world  complexities, 
arriving  at  a  useful  model  for  a  2  terminal  system.   Chapter  I  gives 
the  background  and  a  broad  overview  of  the  communications  network  with- 
in an  Army  division.   The  basic  model,  together  with  some  variations 
are  treated  in  Chapter  II.   More  complex  models,  arrived  at  by  consid- 
ering the  work  of  operators,  are  introduced  in  Chapter  III.   Approx- 
imate solutions  for  the  2  terminal  system,  interconnected  with  a 
k-channel  link,  are  given  in  Chapter  IV.   Chapter  V  provides  an  intui- 
tive insight  into  3  terminal  network  problems.   Areas  for  further 
study  and  extension  of  this  paper  are  discussed  in  Chapter  VI. 


3.   The  Army  Division  Commuaicafcxoris  System. 

There  is  no  rigid  structure  for  the  division  communication  system. 
Many  factors  must  be  weighed  and  evaluated  before  a  concrete  plan  for  a 
communications  system  can  be  drawn  up.   The  type  of  division  (armor, 
airborne,  mechanized,  infantry),  the  mission  of  the  division,  the 
geographical  location,  the  presence  or  absence  of  other  forces,  the 
personnel  strength;  all  can  influence  its  configuration  and  design. 
These  factors  notwithstanding,  there  are  some  general  remarks  which  can 
be  made  about  the  communications  support  given  to  the  division  in  the 
field.   The  division  can  have  from  nine  to  fifteen  battalions.   Under 
present  doctrine,  the  mixture  of  armor,  mechanized  infantry  or  infantry 
battalions  is  tailored  to  the  mission  and  location  of  the  division. 
Control  is  exercised  by  the  Commanding  General  through  three  brigades 
whose  strength  (three  to  five  battalions)  also  is  regulated  by  local 
requirements.   One  brigade  is  normally  held  in  reserve. 

The  signal  battalion,  which  is  organic  to  the  division,  must 
furnish  communication  support  for  the  division  to  Corps  and  Army  HQS 
and  to  the  main  sub-ordinate  headquarters.   Also,  within  the  division 
area  of  responsibility,  signal  centers  must  be  set  up  to  allow  small 
units  to  connect  into  a  division-wide  communications  network.   Many 
modes  of  communication  are  employed.   Hand-held  radios  suffice  at  the 
squad  level;  however  as  the  size  of  the  unit  and  the  physical  separa- 
tion increases,  more  sophisticated  modes  are  required.   These  include 
line  of  sight  microwave  and  tropospheric  scatter.   While  many  advances 
have  been  made  in  the  past  ten  years,  it  is  still  true  that  within  the 
division  the  main  communications  system  is  a  twelve  or  twenty-four 
channel  microwave  radio  system  manually  controlled  by  operators  at 
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switchboards.   Figure  1  shows  how  a  typical  division  system  might  be 
configured.   Usually  every  terminal  has  access  to  Division  Main  and  one 
other  terminal.   This  is  a  precaution  against  loss  of  a  link  due  to 
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12  CHANNEL  RADIO  LINK 


FIGURE  1 


equipment  failure  or  enemy  action.   Division  Main,  normally  the  location 
of  the  Commanding  General  and  his  staff  is  connected  with  as  many  other 
terminals  (at  least  the  brigades  and  division  artillery)  as  equipment 
and  personnel  allow. 

At  each  terminal  an  operator  is  on  duty  to  receive  requests  for 
service,  coordinate  with  other  operators  in  routing  the  call,  and  dis- 
connect the  line  when  the  call  is  complete.   As  shown  in  Figure  2,  each 


TERMINAL 


FIGURE  2 

subscriber  has  an  individual  connection  to  the  switchboard  on  the 
"subscriber"  side  of  the  board.   At  the  switchboard,  a  light  appears 
above  the  termination  of  a  subscribers  line  when  he  turns  a  crank  on 
his  telephone  set.   When  the  operator  responds  to  the  request  for 
service  he  determines  with  whom  the  subscriber  wishes  to  talk.   He  then 
rings  the  operator  at  the  terminal  serving  the  requested  party.   This 
requires  that  at  least  one  channel  in  the  radio  link  between  operators 
is  not  being  used.   When  the  distant  operator  has  reached  the  requested 
party,  the  first  operator  tells  his  subscriber  that  the  connection  has 
been  completed  and  he  can  commence  his  call. 
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CHAPTER  II 

ANALYSIS  OF  A  SIMPLE  NETWORK 
DISREGARDING  OPERATOR  PATCH  TIME 

1.   Basic  Model. 

Before  analyzing  more  complex  models  an  understanding  of  the 
operation  of  a  basic  system  is  helpful.   Accordingly  the  most  trivial 
network  will  be  examined  first.   A  simple  communications  system  is 
shown  in  Figure  3.   Here  we  have  only  one  channel  available  for  use 


1_CHANNEL_ 
RADIO  LINK' 


FIGURE  3 


between  terminal  A  and  B.   The  following  assumptions  will  be  made: 

a)  Requests  for  service  from  the  subscriber  arrive  at  each 
switchboard  with  a  Poisson  distribution,  mean  rate  a. 

b)  The  time  taken  to  complete  a  conversation  from  the  end  of 
routing  to  hang-up,  the  "holding"  time  for  the  call,  is 
distributed  exponentially  with  mean  a • 

c)  Time  between  operator  answering  subscriber  and  the  call 

commencing  (the  operator  "patching"  time)  is  assumed  negli- 

.,.1   •        •         1 
gible  xn  comparison  to  . 

c 

d)  Having  signalled  the  operator,  subscribers  wait  until  the 

operator  answers.   In  other  words,  once  a  subscriber  decides 
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to  make  a  call,  he  is  not  discouraged  by  the  time  required 

to  have  the  operator  answer. 
The  first  two  assumptions  are  reasonable  and  have  been  verified  by  a 
study  on  arrival  rates  and  holding  times  in  a  civilian  telephone 
exchange  in  New  Jersey. 

This  system  can  be  treated  as  a  basic  queueing  model  with  expo- 
nential inter-arrival  time  and  one  server.   The  literature  generally 
denotes  this  as  an  M/M/l  queue.   A  steady-state  solution  for  blocking 
probability  can  be  obtained  if  the  utilization  rate,  p  ,  is  less  than 
unity  where 

p  =  "T"  (1) 

Steady-state  implies  that  the  system  is  in  operation  for  a  long  period 
of  time  and  consequently  the  probability  of  being  in  a  particular  state 
of  the  stochastic  process  is  independent  of  the  state  of  the  process 
when  first  observed.   This  guarantees  that  expected  system  size 
(i.e.,  expected  number  waiting  for  service  plus  those  being  serviced) 
is  always  finite.   The  probability  of  blocking,  p   ,  in  this  system 
is  the  probability  of  the  channel  being  busy  and  is  equal  to 

Pb  "  p  "  T~  (2) 

The  probability  p.   of   (i  -  1)  callers  waiting  for  service  is 

p.   =   (1  -  P)  P1  (3) 

It  is  worth  noting  that  in  this  system,  queue  discipline  is  not 
necessarily  first  come,  first  served  (FCFS) .   Normally  when  queue 
discipline  is  not  specified  it  is  FCFS.   Here,  service  is  not  neces- 
sarily rendered  according  to  order  of  arrival.   For  example,  the  first 
arrival  at  switchboard  A,  with  the  channel  busy,  may  occur  before  the 
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first  arrival  at  switchboard  B.   But  if  the  operator  at  B  becomes 
aware  of  the  termination  of  a  call  before  his  opposite  number  at  A 
then  the  first  subscriber  in  line  at  B  will  be  served  next  although 
some  have  been  waiting  longer  at  A.   This  however  is  an  internal 
organization  problem  and,  while  it  is  of  interest  to  the  individual 
subscriber,  it  does  not  alter  the  blocking  probability. 


2.   Two  terminals  connected  by  a  k-channel  link. 

Consider  now  a  system  consisting  of  two  terminals  with  k  channels 
available  for  communication  between  them.   Figure  4  illustrates  such 
a  system.   The  assumptions  made  in  the  previous  section  will  be 


k   CHANNEL 
RADIO  LINK 


FIGURE  4 


retained.   Implicit  in  these  assumptions  is  the  understanding  that  the 
number  of  lines  on  the  subscriber  side  of  the  board  is  so  large  that 
a)  the  future  arrival  rate,  a  ,  of  calls  is  not  affected  by  the  number 
of  calls  presently  initiated  and  b)  the  probability  of  a  called  party 
being  busy  can  be  assumed  to  be  zero. 

The  system,  therefore,  can  be  modeled  as  an  M/M/k  queue.   That 
is,  inter-arrival  time  exponentially  distributed  with  mean  , 


service  time  exponentially  distributed  with  mean 


,  and  k  servers. 


Define  p   as  the  probability  that  system  size  is  n.   Again  note  that 
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system  size  includes  both  subscribers  waiting  for  service  and  being 
served.   The  probability  p  that  system  size  is  n,  then  can  be  expressed 
as 


n 


( 

2  a 

Po 

a 

) 

n    ! 

n 

2a 

Po 

( 

a 

) 

k  !  k 


n  -  k 


n  <  k 


(A) 


n  >  k 


where   p 


1  -  E  p 
n»l 


The  probability  of  blocking,  p.,  is 


k-1 

1  -  I       p 
n=0 


n 


(5) 


for  this  system. 


3.  Three  terminals  interconnected  by  1  channel  links. 

A  representation  of  this  system  is  shown  in  Figure  5.  Again  the 


1  CHANNEL 
RADIO  LINK 

FIGURE  5 
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assumptions  previously  noted  will  be  retained.   One  additional  require- 
ment for  this  system  is  that  a  call  can  be  routed  over  one  link  only. 
Here  we  denote  the  mean  arrival  rate  at  terminal  i  of  calls  going  to 
terminal  j  as  a..  .   With  operator  patch  time  assumed  zero  and  maximum 
route  length  of  one  link,  the  system  may  be  modeled  by  three  indepen- 
dent M/M/l  queues  with 

a-  +  *-■  (6) 


(7) 


pl 

AB 

BA 

a 

p2 

= 

aAC 

+ 

aCA 

a 

p3 

= 

aBC 

+ 

aCB 

(8) 


As  before  the  order  of  arrival  at  the  switchboard  for  any  of  these 
queues  is  not  necessarily  the  order  of  service. 

4.   Three  terminal  system  with  alternate  routing  available. 

Introduction  of  alternate  routing  gives  a  degree  of  flexibility 

to  the  previous  model  as  well  as  adding  realism.   Alternate  routing 

allows  a  call  from  A  to  B  to  be  routed  on  links  A-C  and  C-B,  if  they 

are  free  and  link  A-B  is  busy.   Here  assumptions  made  in  Section  1 

above  hold  with  the  exception  that  if  all  links  are  busy  incoming 

callers  hang  up.   For  ease  in  computation  assume  all  arrival  rates, 

a   ,   a^.,   .  .  .  ,   a^   are  the  same  and  equal  to   a/2.   Figure  6 
AB     BA  LB 

depicts  all  possible  states  of  the  system.   A  state  is  a  unique 
position  which  completely  describes  the  process  at  any  instant  of  time 

It  is  now  possible  to  obtain  a  transition  rate  matrix  A,  which 
gives  the  rate  of  movement  between  states  of  the  system.   Solving  the 
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all  lines  empty 


one  line  busy 
1  call 


& 


S *"  c 


2  lines  busy 
2  calls 


2  lines  busy 
1  call 


3  lines  busy 
2  calls 


3  lines  busy, 
3  calls 


POSSIBLE  STATES  OF  THE  SYSTEM 


FIGURE  6 
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system  of  equations 

P  A   =   0  (9) 

->- 
where  P  =    (p  p   .  .  .   p  ,)    p.,  the  steady  state  probability 

that  the  system  is  in  state   i  can  be  obtained.   Figure  7  shows  the 

transition  matrix  for  this  system.   Note  that  it  is  possible  to  reduce 

the  transition  matrix  from  14  x  14   to   6  x  6  by  grouping  similar 

states  together  in  one  "super-state". 

For  example:   states  2,  3,  and  4  may  be  combined  to  create  the 

super-state  —  "one  call,  one  line  busy".   Using  equation  (8),  a  set 

of  simultaneous  equations  can  be  written  as  follows: 


■a  px     +  o   pn  +  o   pIV  =  0 

a.  Pl     -  (a  +  a)   P][I     +       -y-  a  pm  +  ~  PV  =  ° 


y-a  pXI     -(a  +  yo)Pm  +  a   p  =   0 


(a  +  3   }    PIV     +"TPV  =  ° 


+     T-     Ptu      "  T  °Pu  =   ° 


(10) 


3        FII  '3        FIV  3   wpV 


a  PIII  "  a  PVI  =   ° 


PI  +  PII  +  PIII     +  PIV       +        PV     +     PVI  =   x 
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FIGURE    7 


TRANSITION      MATRIX 
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a       1 
As  an  example,  let  =  — ~—  .   This  yields  the  following  results 


p     =  probability  of  no  lines  busy 


112 


I       v 3  "J  221 


p     =  probability  of   1  line  busy 


52 


n      r*™— j  —   -  ~™  —.3  221 


(11) 


p    +  p     =  probability  of   2  lines  busy 


p..   +  p     =  probability  of   3  lines  busy 


30 


221 


27 


v    ,  Fvi      r~— -— j  w*   -  ..—  — ,  2n 


a      1 
For  this  system  the  probability  of  blocking,  p,  ,  with  =  — —  is 

2 
p,   =  prob   {3  lines  busy}  +  — —     prob   {2  lines  busy} 

27        20  57  ^12^ 

+  -SBHP-   =   -^4t—     =    .259 


221       221         221 

For  the  same  system  with  2  channels  per  link  connecting  every  pair  of 
terminals,  there  are  34  possible  states.   Using  symmetry  and  grouping 
similar  states  as  before,  a  solution  can  be  obtained  from  a  set  of 
11  simultaneous  equations.   It  can  be  recognized  that  even  this  simple 
problem  requires  a  significant  amount  of  time  to  solve  manually  if 
computing  equipment  is  not  available. 

5 .   Extensions  -to  multi-terminal  networks 

It  would  be  useful  if  generalizations  of  the  results  of  the  previous 
section  were  available  for  n-terminal  networks  interconnected  with 
k-channel  links.   No  such  solution  seems  to  be  available.   It  appears 
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that  systems  of  interconnected  terminals  must  be  treated  individually 

and,  except  for  the  most  trivial,  the  combinatorial  aspects  make 

[3] 
analytic  solutions  appear  unpalatable.   Benes      states  "the  perfor- 
mance of  a  connecting  network  for  a  given  level  of  offered  traffic  is 
determined  largely  by  its  configuration  or  structure  ...".   This 
statement  certainly  has  intuitive  appeal  and  an  example  can  exhibit 
its  validity.   Consider  two  networks,  with  the  same  number  of  terminals 
and  total  channels  available  for  use,  as  shown  in  Figure  8.   They 
differ  only  in  network  configuration.   Nevertheless,  certain  states 


2  \1 

1  CHANNEL 


D 


RADIO  LINK 
(ii) 


FIGURE  8 

will  have  different  probabilities.   For  example  consider  the  state 
"2  lines  busy,  1  call  ".   For  (i)  there  are  eight  ways  to  obtain  this 
condition  while  in  (ii)  there  are  only  six  ways  to  accomplish  the  same 
condition  even  though  the  number  of  terminals  and  total  available 
channels  are  the  same.   Obviously  with  different  flow  rates  the  systems 
would  be  incompatible. 
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CHAPTER   III 


MODELS  FOR  NETWORKS  WITH  NON-ZERO  PATCH  TIME 


1.   Two  terminal  system  with  operator  patch  time. 

A  more  complex  system  to  analyze  is  one  similar  to  that  shown  in 
Figure  3  but  in  which  operator  patch  time  is  considered  significant. 
Again,  we  assume  mean  arrival  rates  of  calls  at  both  switchboards 
of  a  and  holding  time  for  calls  with  mean  rate  o.   In  addition  operator 
patch  time  is  taken  as  distributed  exponentially  with  parameter  y. 
Queue  discipline  is  nominally  first  come,  first  serve  and  callers  wait 
for  service  if  the  channel  is  occupied. 

Having  only  a  one  channel  link  between  terminals  makes  available 
for  this  system  a  special  analytical  solution.   Service  time  is  the 
sum  of  two  independent  random  variables  distributed  exponentially  but 
with  different  means.   Normally  once  an  operator  has  completed  a  patch 
he  is  free  to  assist  the  next  caller  in  line  at  the  switchboard.   In 
this  system,  however,  as  there  is  only  one  channel  the  operator  is 
idle  until  the  call  occupying  the  channel  is  terminated.   In  practice 
the  operators  arrange  the  routing  of  the  requested  call  on  the  same 
channel  that  the  call  itself  will  use. 

The  system  can  be  modeled  as  an  M/G/l  queue  with  arrival  rate 
distributed  exponentially  with  mean  2a.  Service  time,  B(t),  has  a 
mean  and  variance 

m   =   -i-  +   JL    =      °  +  ^  (1) 

o  \i  a\i 

v  =   -L-  +   -i-    =    q2  +  ^  (2) 

2  2  2  2 
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[4] 
The  embedded  Markov  chain  technique  of  Kendall    can  be  used  to  solve 

this  type  of  queue  and  obtain  answers  to  the  most  interesting  questions 

(e.g.   number  of  callers  awaiting  service,  probability  of  having  x 

callers  awaiting  service,  expected  waiting  time  for  service) .   Define 

X   ,   n   0  the  number  of  calls  awaiting  service  when 

n  ntn  call  is  completed 

U   ,   n  >  0  the  number  of  calls  arriving  at  the 

n 

switchboards  during  the  routing  and 

calling  of  the  n  +  1st  call. 

u.  ,   k  >    0  prob  {U  =  k} 

k      =  n 


p  traffic  intensity 

5 

y 


G   (z)  z  transform  of  y 


F   (s)  Laplace  transform  of  y 

y.  steady  state  probability  that  x  =  i 

We  know     that   {x  }  is  a  recurrent  Markov  chain  if 

n 

G  *  (1)   -  Z   k  u.    <,  1  (3) 

Uk  k=0     k 


but, 

i 

G  *  (1)    =   E[U  ]    =  a  m   =  p  (4) 

k 

and  so  if  p  is  less  than  unity,  steady-state  solutions  exist.   It 

can  be  shown  that 

(1  -  p)  (1  -  z)   F   *   [2a  (1  -  z)] 
G    (z)    =  (5) 

1  FB(t)  [2a  (1  "  Z)]   "  Z 


where 


FB(t)  [2a  (1  -  z)]   =  (6) 


[a  +  2a  (1  -  z)]  [y   +  2a  (1  -  z)] 

in  this  example,  and 
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yQ   =   G     (z)/z=0   =    1  -  p 
i 

y1     -   G  *  (z)/z=0  (7) 


*x 


y.   =   G    (z)/z=0  ,  etc. 
i       y . 
l 


To  illustrate  the  solution  technique  the  following  example  is 
given : 

let   a   =   — s—    per  hour 

y   =   50     per  hour 

o   =   10     per  hour 

.           _  .   50  +  10         6        . 
then  p   =   5   5QQ =   -j^—       <   1 

and  steady  state  conditions  exist. 

y   =   1  -   p   =   .4 
J  o  K 

p,   =   1  -  y      .6 

rb  J  o 

Prob  {  more  than  1  call  at  the  switchboard)  =  1  -  y  -  y1 
=   1  -  .4  -  .23   =   .37 

Expected  system  size  is  obtained  from 

E[x  ]    =    P  +  P2  +        "2  v 
n  2(1  -  p) 

'  ^  "   1'375 
Operator  idle  time  per  hour     y   H r—  y1 

=   .515   hours 
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Expected  total  time  in  system   =  ■= —  E(x  ) 

2a     n 

=   .275   hours 

Expected  time  waiting  for  an  operator  =   .275  -  m 

=   .156  hours 

2.   Two  terminal  system  with  one  way  traffic  flow. 

A  much  more  difficult  system  to  analyze  is  that  shown  in  Figure  4 
with  the  following  additional  constraints.   All  traffic  is  assumed  to 
be  initiated  at  one  terminal  and  exponential  operator  patch  time  (with 
mean  rate  u  )  is  considered.   Now  once  the  operator  has  completed  a 
patch  from  A  to  B  he  can  immediately  begin  to  service  the  next  customer 
(if  there  is  one)  at  the  switchboard  providing  at  least  one  of  the  k 
lines  is  vacant. 

This  system  can  be  described  as  a  two  stage  queue.  Stage  1  is  the 
operator  patching  operation  and  is  an  M/M/l  queue  with  infinite  waiting 
room.  Stage  2  is  the  holding  time  per  call  and  is  a  G/M/n,  n<  k,  queue 
with  no  waiting  permitted. 

The  literature  includes  work  done  in  the  area  of  sequential  array 
or  tandem  queues.    '     This  system,  however,  is  unlike  those  treated 
in  the  literature  in  that  the  stages  are  not  alike.   They  have  a  strong 
dependence  on  each  other,  unlike  the  usual  tandem  queues  which  involve 
movement   from  one  stage  to  the  next  and  where  the  only  dependence 
of  "successor"  stages  is  upon  their  predecessor's  output.   This  system 
has  the  feature  that  while  Stage  2  takes  its  arrival  rate  as  the 
service  output  of  Stage  1,  the  service  rate  in  Stage  1  also  depends 
upon  the  system  size  of  Stage  2. 
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A  solution  technique  for  this  system  is  as  follows: 
Let  r  denote  the  number  of  channels  busy  (in  either  patching  or 
conversation).   Then  for  Stage  1,  service  rate  for  a  customer  is 


Hunt     shows  that  for  a  2  stage  system,  both  stages  being  M/M/l 
with  unlimited  waiting  room,  the  probability  of  system  size  being  n  is 

P  (z)   -  I  pa)  p  ^ 

*n        .  n    i   rn-i 
i=0 


=   (n  +  1)  (1  -  p)2  pn 
where  p.     =  probability  {i  calls  in  stage  j  at  steady  state} 
p.(R)   =  probability  {i  calls  in  an  R-stage  system} 

Stage  1  is  recurrent  (i.e.,   {p.   }  ,  _i  >_ 0   is  a  recurrent  Markov 
chain)   if  p,  the  traffic  intensity  or  utilization  rate,  is  strictly 

less  than  unity.   Normally  in  an  M/M/l  queue  p  =  ,  however,  in 

this  case  service  rate  depends  upon  whether  any  lines  are  empty.   We 
can  therefore  approximate  u,  the  expected  mean  service  rate  of  Stage  1 
to  be 

u  =  u*  prob  {r  <  k}   +  0-  prob  {r=k}  (9) 

=  u-  prob  {r  <  k} 
a  constant. 


Stage  1  is  recurrent  if  a  <u.   Then  the  following  statement  is  true: 

If  Stage  1  is  a  recurrent  Markov  chain  the  rate  of  departure  from 

T8 1 
Stage  1  must  equal  the  arrival  rate  a  .      The  distribution  of  the 
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inter-departure  times  is  no  longer  exponential  due  to  the  change  in 
service  time  distribution.   However,  for  many  problems  of  interest  the 
output  of  Stage  1  can  be  approximated  as  Poisson  without  incurring 
great  error.   This  is  because  in  any  real-life  system  the  prob  j  //  of 
lines  busy  =  k  j-  is  small. 

A  series  of  simulations  were  run  on  an  IBM  360/67  computer  to 

[9] 
verify  this  Poisson  approximation.   An  internal  package  called  GPSS 

was  used  to  generate  3  runs  of  a  2  stage  system  similar  to  that 

described  above.   1000   inputs  were  used  in  each  run  and  reasonable 

values  of  a,  y  and  a   were  selected.   Data  from  these  simulations  are 

shown  in  Table  1.   It  can  be  seen  that  the  inter  departure  times  of 

the  output  of  Stage  1,  a   have  the  general  characteristics  of  an 

exponential  distribution  (i.e.,  mean  and  standard  deviation  are  the 

same) .   In  addition  the  mean  rates  of  input  and  output  of  Stage  1  are 

approximately  the  same.   By  way  of  contrast,  2  runs  were  generated  with 

parameter  values  selected  so  that  the  system  was  unstable.   Runs  #4 

and  //5  of  Table  1  show  the  difference  in  input  and  output  ratio  in 

Stage  1  when  the  service  rate  of  Stage  2  is  equal  to  input  rate  of 

Stage  1  (run  #4)  and  when  the  service  rate  of  Stage  2  is  less  than  the 

arrival  rate  to  Stage  1  (run  #5) . 

In  Stage  2  the  probability  of  having  2  in  the  system  is 

(2)        1  p2± 

pu;  =  ~ — —  (10) 


i 


k        j    i! 

I      p2 


j=0    -, 


l 

QJ_ 

v     J 
i!  .En  p2 

J=0  IT 
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PARAMETERS 


RVN 

111 
a,          u            a 

n 

0 

mean 

■2 

stand. 

dev. 

1 

10 

2 

10 

6 

9.74 

9.79 

2 

3 

2 

10 

6 

2.61 

2.69 

..  .  _, 

^ 

3 

5 

2 

10 

6 

4.45 |    4.56 

J , 

4 

5 

1 

30 

6 

T     1 

5.94     5.88 

j 

5 

1     3 

1 

40      { 

6 

6.84     7.02  f 

i 

I 

1                           1 

arrival  rate 


arrival  rate  from  stage  1 


ervice  rate 


service  rate 


STAGE  1 


STAGE  2 


NOTE:   numbers  in  table  represent  time  units  between 
arrivals,  departures  or  services 


RESULTS  OF  SIMULATION  OF  2-STAGE  QUEUE 


TABLE  1 
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k 
and  Prob   {r=k}   =    P2 (11) 

i  i   k       J 
k!   z     p2j 


J-0 


j! 


where   pn  =  — - 
2    c 


In  Stage  1  the  system  utilization  rate  is 


P 


1  k 

u  -  u  p0 


k! 

j-0       3' 


a   k!       E  P2 


3=0        j! 


k 

u  k!       E      p„J      -  up„ 
j=0— l— 


n  k 


i=0 
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(12) 


and  p<1}      =      (1   -   P]_)    p1i  (13) 


Therefore  for  the  2  Stage  system 


n    D(1)   D    (2) 

E   Pi   Pn- i  n  <  k 

i=0 


and  by  using  equations  8  through  14  we  get 

n-i 
n  .  p^ 

i£o    &  ~  Pi>  pi1   k 53    ,n  *  k 


(n-i)    ! 


j=0       j! 

,<»  - 

k  k-i 

n         s         (n-k+i)  p2 

i£o    (1-pi}  pi  k Fp-  'n  =  k 

(k  -  i)    !j=0     ~7T 


and  .  i        k-i 

K  P-i  Po 

Pb     -      (1  "  pi>    i£o    k pTJ-  (16) 

x.       2 


(k-i)!      j=0        j 


is  the  probability  of  blocking, 
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CHAPTER   IV 
2  TERMINAL  SYSTEMS  WITH  PRIORITIES  AND  UNLIMITED  ACCESS 

1.   Models  with  unlimited  access. 

A  much  more  complex  system  to  analyze  is  that  in  Figure  A  when 
arrivals  from  subscribers  are  allowed  at  both  switchboards.   As 
discussed  in  Section  2,  Chapter  III  the  system  can  be  modeled  as  a 
tandem  network  consisting  of  a  M/M/l  queue  and  a  G/M/k  queue  with  no 
waiting  room.   A  subtle  distinction  should  now  be  noted.   The  model  is 
now  susceptible  to  total  collapse  under  the  following  circumstances. 
Having  assumed  Poisson  distributed  arrivals  from  subscribers  at  both 
switchboards,  the  probability  of  instantaneous  arrivals  at  both  switch- 
boards is  zero.   However,  consider  what  would  happen  if  the  time 
between  arrivals  is  small  compared  to  the  operator  patch  time.   Operator 
B  could  have  begun  patching  a  subscriber  before  operator  A  could  patch 
his  subscriber  to  B.   Thus  the  system  would  be  suspended  in  a  state 
of  inaction  with  both  operators  trying  to  contact  each  other.   Some 
real  world  networks  prevent  this  situation  from  occuring  by  designating 
one  channel  for  use  by  operators.   Hence  both  operators  would  never  be 
trying  to  get  in  touch  with  each  other  at  the  same  time.   Use  of  the 
"orderwire",  as  the  designated  line  is  sometimes  called,  is  inefficient 
because  only  k-1  channels  can  be  used  for  communication  from  A  to  B. 
In  most  systems,  if  the  problem  arises,  it  is  disposed  of  by  putting  a 
higher  priority  on  calls  arriving  at  the  switchboard  from  another 
operator  than  on  calls  from  the  subscriber  side  of  the  switchboard. 
If  an  operator  is  busy  with  one  of  his  subscribers  and  notices  the 
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other  operator  ringing  he  immediately  stops  service  and  attends  to  the 
priority  call. 

Since  the  mean  rate  of  calls  at  both  switchboards  is  assumed  to 

and 


be  — j-   ,  each  switchboard  in  effect  receives  — r—  high  priority  , 
low  priority  calling  rates.   Operator  patching  time  can  no  longer  be 
considered  an  integral  block  but  must  be  reduced  to  its  two  component 
phases:   the  work  done  by  each  operator-   Without  much  error,  their 

u) .   The  system  can  then 
be  depicted  as  shown  in  Figure  9.   Note  that  exit  from  queues  A  and  B 


rates  can  be  considered  equal  (uA   =  u„ 

n      A      B 


SWITCHBOARD 
A 


G/M/k 

a 

i 

L 

5 

a 

PRIORITY 


PRIORITY 


G/M/l 
u 


SWITCHBOARD 
B 


FIGURE   9 

at   1   and   2   respectively  is  allowed  only  if  the  number  being 
served  at   3   ,    4    and    5     plus  the  number  waiting  at   3 
and   4   is  less  than  k.   Queue  C  has  no  waiting (as  indicated  by  k 
servers) .   Once  entry  is  obtained  to  queue  C   service  begins  immedi- 
ately.  This  corresponds  to  the  fact  that  a  call  is  not  routed  unless 
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an  empty  channel  is  available.   Once  the  routing  is  complete  the  call 
can  begin  (queue  C)  without  delay. 

2.   An  approximation  of  the  system. 

Analysis  of  this  system  is  extremely  complex.   It  combines  in 
varying  degree  such  features  as  tandem  queues  with  restricted  waiting 
facilities,  networks  of  queues,  and  priority  queues  with  general 
arrival  distributions.   All  of  these  topics  have  been  discussed  individ- 
ually in  the  literature  and  solutions  for  some  of  the  more  useful 
properties  obtained. 

The  following  approximation  would  be  helpful  in  obtaining  blocking 
probabilities  for  this  system:   Assume  that  the  inter  departure  time 
from  the  queues  for  patching  is  distributed  exponentially.   It  is 
obvious  that  in  steady  state  the  output  of  the  first  queue  (operator) 
is  no  longer  Poisson.   Although  arrivals  from  the  subscriber  side  enter 
at  a  Poisson  rate,  because  they  have  lower  priority  and  because  service 
is  sometimes  interrupted  when  k  lines  are  full,  they  leave  at  a  mean 
rate  of   9   but  with  unknown,  general  distribution.   If,  however,  the 
total  arrival  rate  a  is  small  compared  to  u  (effective  service  rate, 
see  Section  2,  Chapter  III)  the  distribution  of  outputs  entering  the 
second  switchboard  queue  should  be  approximately  Poisson. 

To  show  that  this  approximation  can  be  very  reasonable,  a  queue 
having  both  low  priority  and  high  priority  arrivals  was  simulated. 
Both  types  of  arrivals  were  Poisson  with  mean  rate  15  per  hour.   Service 
rate  was  exponential  with  mean  time  of  1  minute.   The  output  distri- 
bution is  shown  in  Figure  10.   Because  of  the  low  utilization  rate,  0.5, 
the  inter  departure  time  of  the  low  priority  calls  has  the  shape  of  an 
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exponential . 

If  this  approximation  is  used  the  system  can  be  modeled  as  in 
Figure  11.   It  is  clear  that  only  one  operator  can  be  busy  handling 
calls  from  the  line  side  of  the  switchboard  at  any  one  time.   Assuming 
that  arrivals  at  the  switchboard  from  the  line  side  are  Poisson 


STAGE  1 


STAGE  2 


FIGURE  11 

distributed,  entry  into  the  waiting  line  for  Stage  1  (call  transfered 
from  one  operator  to  another)  is  permitted  only  if  system  size  (number 
being  served  in  Stage  2  plus  number  waiting  in  Stage  1  plus  1  or  0 
being  served  in  Stage  1)  is  less  than  k. 

3.   System  State  Diagram. 

Define  p..   to  be  the  steady-state  probability  of  having  i  calls 
waiting  to  be  patched  or  being  patched  by  the  operator  and  j  calls  in 
progress.   The  state  diagram  for  this  system  is  shown  in  Figure  12 
(note  that  i  +  j  cannot  be  greater  than  k) .   The  probability  of  r 
lines  being  busy  (where  0  <_  r  <_  k)  is  the  sum  of  the  state  proba- 
bilities of  the  diagonal  where  i  +  j  =  r. 
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4.   System  Transition  Matrix. 

A  precise,  general  analytical  solution  of  this  system  does  not 
appear  possible.   This  will  be  more  apparent  later.   The  state  diagram 
can  be  broken  into  k  +  1  "super-states",  each  super-state  being  the 
probability  of  r  lines,  0  <_  r  <.  k,  being  busy.   While  all  states 
within  a  super-state  have  the  same  probability,  aAt,of  moving  to  the 
next  higher  super-state  in  At,  each  has  a  different  rate  of  movement 
to  a  lower  state.   Transition  within  states  of  a  super-state  is  at 
equal  rates  u.   Thus  the  probability  of  going  from  super-state  r  to 
r  -  1  is  dependent  on  which  particular  state  of  the  super-state  is  the 
departure  point.   A  transition  matrix  for  rate  of  movement  between 
super-states  can  be  obtained  as  shown  in  Figure  13.   The  d   terms  are 
unknown  constants,  which  could  possibly  depend  upon  a, a  and  u. 

The  number  of  lines  occupied  is  a  "birth  and  death  process"    ; 
movement  is  permitted  only  to  adjacent  states  of  the  stochastic 
process  at  any  instant  of  time.   An  expression  for  tt  ,  the  steady- 
state  probability  of  having  r  lines  busy  regardless  of  the  state  of  the 
system  when  first  observed  is 


1  <  r  <.  k       (1) 
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5.   Determination  of  f  . 

r 


The  matrix  above  was  solved  by  computer  for  6  sets  of  parameters. 

-*■ 

From  the  results  obtained  for  the  vector  tt  ,  values  for  f  were 

r 

obtained  by  noting  (2)  can  be  written 


£(+) 


fr   =  —?-      (-—)  (3) 


Figure  14  shows  a  plot  of  f   against  r  for  these  solutions. 

6.   Error  introduced  by  approximations. 

To  test  the  solutions  obtained  with  equation  4-2  by  using  the 
approximations  mentioned  previously  (i.e.,  assuming  blocking  does 
not  seriously  affect  the  distribution  of  inter  departure  times  and 
disregarding  the  effects  of  priority  arrivals),  systems  were  simulated 
using  the  same  parameters.   Figures  15  through  17  compare  the  exact 
(simulated)  probabilities  of  channels  being  used  with  the  approximate 
solutions.   In  addition,  a  sensitivity  analysis  was  performed  on  the 
approximate  solutions  to  test  the  variation  in  blocking  probability 
with  changes  in  a.   These  results  are  shown  in  Figure  18.   The  impli- 
cations of  these  results  are  far-reaching.   It  appears  that  the 
probability  of  blocking  is  a  linear  function  of  arrival  rate,  a. 
The  upper  boundary  on  a  would  be  the  minimum  of  patching  rate  or  k 
times  the  calling  rate.   This  information  would  be  useful  to  the 
systems  planner  who  is  confronted  with  fluctuations  in  arrival  rate 
over  time. 
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CHAPTER  V 
3  TERMINAL  NETWORKS 

1.  Simulation  results. 

A  simulation  of  a  3-terminal  network  was  run  to  determine  the 
effects  on  blocking  probability  of  allowing  routing  from  terminal  x 
through  terminal  y  to  terminal  z  when  all  k  lines  between  x  and  z  are 
busy.   Intuition  would  indicate  that  the  mean  number  of  channels  busy 
for  each  pair  of  the  3-terminal  network  would  be  higher  than  that  for 
a  2-terminal  network  with  the  same  parameters.   This  is  due  to  the 
inefficiency  of  routing  (i.e.,  using  2  channels  in  tandem  to  route  one 
call  when  the  direct  route  is  blocked)  in  the  3-terminal  system. 
However,  it  could  be  expected  that  the  probability  of  blocking  for  any 
call  would  be  substantially  reduced  because  at  least  2  links  must  be 

at  full  capacity  to  block  a  call.   This  would  happen  despite  higher 

-> 
individual  link  blocking  probabilities.   Table  2  lists  the  tt  vector 

for  each  link  of  the  3  terminal  system  and  that  of  a  2  terminal  system 

with  the  same  parameters.   Results  from  the  simulation  showed  that 

78%  of  calls  in  3-terminal  networks  which  were  initially  blocked  could 

be  indirectly  routed  to  their  destination  immediately. 

2.  n  -  terminal  networks  with  alternate  routing. 

It  would  appear  that  an  upper  limit  exists  for  the  amount  of 
alternate  routing  that  can  be  used  in  a  general,  n-terminal  network. 
The  inefficiency  of  the  system  should  increase,  as  a  function  of  n, 


44 


2-TERMINAL    3-TERMINAL    3-TERMINAL    3-TERMINAL 
SYSTEM        LINK  1        LINK  2        LINK  3 


:0 


'10 


11 


12 


.0019 
.0049 
.0319 
.0618 
.0917 
.1565 
.1525 
.1326 
.1276 
.0737 
.0568 
.0568 
.0508 


.0014 
.0128 
.0381 
.0648 
.1109 
.1406 
.1401 
.1396 
.1194 
.0965 
.0584 
.0595 
.0521 


.0009 
.0099 
.0298 
.0690 
.1202 
.1361 
.1376 
.1411 
.1182 
.0725 
.0621 
.0511 
.0506 


.0030 
.0100 
.0356 
.0658 
.1130 
.1391 
.1386 
.1396 
.1140 
.0884 
.0618 
.0591 
.0512 


TABLE   2 
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until  the  probability  of  blocking 


pb  =    v  ■  v  •  v  •  •  '  \   ,         (1) 

12     3         n-1 


(where  tt,    is  the  probability  of  k  channels  busy  in  link  1) 
kl 

is  greater  than  the  probability  of  blocking  for  a  2  terminal  network 
with  the  same  parameters.   Efficiency  is  defined  as  the  ability  of  the 
network  to  route  calls  over  the  minimum  number  of  links. 


4b 


CHAPTER  VI 


ANALYSIS   OF  RESULTS 


1.   General. 


A  study  of  Figures  15,  16,  17  indicates  that  the  approximate 
solution  for  the  distribution  of  number  of  lines  busy  in  a  2^terminal 
system  is  satisfactory  for  many  purposes.   While  the  absolute  differ- 
ence between  the  actual  and  approximate  results  is  small,  the  differ- 
ences are  significant  and  the  approximate-  solutions  should  be  used 
for  systems  design  or  operational  planning  only  if  simulation  methods 
are  unavailable.   The  approximate  results  serve  a  purpose;  they  provide 
"ball-park"  answers  for  the  systems  planner  when  formerly  no  infor- 
mation appears  to  have  been  available. 

The  "real  world"  condition  that  after  k  lines  are  busy  each 
operator  can  service  one  more  subscriber  was  neglected  in  the  approx- 
imate solution.   This  last  service  at  the  switchboard  remains  there 
until  one  of  the  k  lines  is  free.   Then  it  is  patched  to  the  other 
operator  and  frees  the  queue  from  which  it  came.   This  would  have 
added  another  diagonal  row  of  states  to  the  state  diagram,  Figure  12. 
Consequently  these  state  probabilities  would  be  added  to  p   and  the 
approximate  solution  should  be  closer  to  the  actual. 

The  function  f   (Figure  14)  appears  to  be  very  stable  for  the 
range  of  a,  a  and  u  for  which  the  approximations  are  valid.   Figure 
14  which  shows  that  f   is  linear  over  the  range  of  r  is  a  very  impor- 
tant discovery.   From  this  curve,  the  probability  of  r  channels  busy 
can  be  determined  by  equation  4-2  if  o  and  a  are  knovJn.   Extrapolation 
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for  larger  values  of  k  would  be  possible.   It  is  interesting  to  note 
the  similarity  between  the  results  of  the  2- terminal,  k-channel 
networks  with  and  without  operator  patch-time  considered.   Equations 
4-2  (with  patch-time)  and  2-3  differ  only  by  constants.   This  is 
reasonable  because  the  patch-time  might  be  considered  as  just  a  portion 
of  a  larger  unit  of  service  time  (patch-time  and  call  holding-time) . 
Consequently  the  systems  differ  only  in  mean  service  time. 

2.   Areas  for  future  investigation. 

Among  some  of  the  more  fertile  areas  for  further  study  would  be 
an  analytical  solution  for  the  2  terminal,  k-channel  network  without 
availing  of  the  approximations  used  in  this  study.   Introduction  of 
priority  arrivals  from  the  subscribers  into  the  models  used  here 
would  also  be  a  useful  extension.   It  is  doubtful  that  networks  with 
more  than  2  terminals  would  yield  solutions  analytically  without 
benefit  of  over-simplifying  assumptions.   The  nature  of  such  assump- 
tions would  also  provide  a  fruitful  investigation. 
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